Dissipative Quantum Systems and the Heat Capacity Enigma 
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We present a detailed study of the quantum dissipative dynamics of a charged particle in a 
magnetic field. Our focus of attention is the effect of dissipation on the low- and high-temperature 
behavior of the specific heat at constant volume. After providing a brief overview of two distinct 
approaches to the statistical mechanics of dissipative quantum systems, viz., the ensemble approach 
of Gibbs and the quantum Brownian motion approach due to Einstein, we present exact analyses 
of the specific heat. While the low-temperature expressions for the specific heat, based on the 
two approaches, are in conformity with power-law temperature-dependence, predicted by the third 
law of thermodynamics, and the high-temperature expressions are in agreement with the classical 
equipartition theorem, there are surprising differences between the dependencies of the specific heat 
on different parameters in the theory, when calculations are done from these two distinct methods. 
In particular, we find puzzling influences of boundary-confinement and the bath-induced spectral 
cutoff frequency. Further, when it comes to the issue of approach to equilibrium, based on the 
Einstein method, the way the asymptotic limit (t — > oo) is taken, seems to assume significance. 

PACS numbers: 05.70. -a, 05.30. -d, 05.40. Jc 



I. INTRODUCTION 



Recent years have seen great strides in the statistical 
mechanics of dissipative quantum systems Dissipa- 
tion arises when the quantum degrees of freedom of a 
heat bath, which is strongly coupled to a subsystem of 
interest, are projected (or integrated) out of the Hilbcrt 
space of the total system. Two different approaches, de- 
tailed below in Sec. II, have been used in this context: (i) 
the usual Gibbs approach that focuses on the partition 
function [2| and (ii) the Einstein approach that hinges 
on a quantum Langevin equation for the subsytem [3(. 
Lately it has been argued that the presence of quantum 
dissipation yields a satisfactory behavior of the funda- 
mental thermodynamic attribute, viz., the heat capacity, 
as far as the low-temperature properties are concerned 
Here we will point out that there are some puzzling 
issues even for the high temperature limit of the heat 
capacity, apart from the intriguing low-temperature at- 
tributes. Before we address this question, it is important 
to review the kind of subsytem we have in mind and the 
foundational basis of statistical mechanics, which we do 
below. While our present discussion as well as that in 
Sec. II are set within the domain of classical statistical 
mechanics, extension to quantum mechanics can be eas- 
ily carried out, as indicated in Sec. III. But we want to 
first concentrate on some preliminaries about the subject 
of statistical mechanics itself. 

Statistical Mechanics provides the microscopic basis of 
the macroscopic poperties of a system described by the 
subject of thermodynamics. Though the power of sta- 
tistical mechanics comes to the fore in its full glory for 
an interacting many body system, such as in the exact 
formulation of second order phase transitions by means 
of the two-dimensional Ising model [f| , many of the in- 
tricacies can be elucidated for just a single entity, albeit 



in contact with a heat bath comprising an infinitely large 
number of (invisible) degrees of freedom. It is this sim- 
plified approach to statistical mechanics in the context of 
a single particle embedded in a heat bath that we shall 
adopt in this paper. 

The dynamics of a particle of mass m is described by 
the system Hamiltonian defined by 



P 2 



(i) 



where p is the canonical momentum vector of the parti- 
cle moving under an arbitrary potential V(q) which is a 
function of the generalized coordinate vector q. We shall 
discuss three distinct cases in the sequel: 



(a) Free particle: 

V{q) = , 

(b) Harmonic oscillator : 

V(q) = imwgf 2 , 
loq being the frequency of the oscillator, and 



(2) 



(3) 



(c)Charged oscillator in a magnetic field, that is de- 
scribed by a momentum and coordinate-dependent po- 
tential: 



V(q,p) = - (p.A(q)+A(q).p) 
Zmc 



(4) 



A(q) being the vector potential, the curl of which yields 
the magnetic field B: 



fi = Vx A(q) 



(5) 



2 



It is evident that for zero vector potential, case (c) re- 
duces to (b). If additionally, uj is also zero, case (a) is 
obtained. In what way are these limiting situations ar- 
rived at, for a quantum dissipative system, will indeed 
be the focus of our discussion below. 

It should be mentioned here that the problem of a 
charged oscillator in a magnetic field is relevant in the 
context of Landau diamagnetism [6] which has had a deep 
impact on modern condensed matter physics through 
phenomena such as the quantum Hall effect Q. Lan- 
dau diamagnetism, which is purely quantum in origin, 
is characterized by strong boundary effects that can be 
mimicked by the oscillator potential Q . The presence of 
a quantum bath, comprising of, say, bosonic excitations 
like phonons, lends additional richness to the problem as 
it allows us to study the effect of dissipation on Landau 
diamagnetism [9(. In this article however our focus of 
attention is not diamagnetism but the thermodynamic 
property of the heat capacity. 

The microstate of the particle at a given time is speci- 
fied by a point in the 6-dimensional (three for coordinates 
and three for momenta) phase space. As the time evolves 
the phase point curves out a phase trajectory. While in 
classical mechanics the trajectory is uniquely determinis- 
tic, once the initial values of q and p are given, the point 
of statistical mechanics is that the phase trajectory ran- 
domly changes from one 'realization' of the system to an- 
other. The meaning of 'realization' becomes clear if one 
considers how experiments are performed. A realization 
corresponds to a given experiment when one watches the 
trajectory evolve in time. Of, course, the whole statisti- 
cal basis of data collection is to repeat the experiment, 
this time tracking a different trajectory, even though the 
initial values of {<f, p] are the same. It is this multitude 
of trajectories corresponding to multiple realizations of 
the system that yields the concept of 'ensemble' in sta- 
tistical mechanics — an ensemble means a collection of 
possible realizations of the system. Thermal equillibrium 
is said to be reached when experiments are repeated so 
many times that all possible trajectories (realizations) in 
the phas e sp ace are explored-this yields the notion of 
'mixing' [10| . 

With these preliminaries the outline of the paper is 
as follows. In Sec. II, we review the Gibbs and Einstein 
approaches to statistical mechanics. Although our treat- 
ments are couched in classical terms similiar results hold 
for quantum phenomena as well. With these approaches 
in the background we summarize in Sec. Ill, the newly 
developed subject of dissipative quantum systems. In 
Sec. IV we analyze the results for the heat capacity for 
the three problems (a-c) and point out certain surprises 
when we consider the various limits of case (c). In Sec.V, 
we summarize the results. 



II. GIBBS AND EINSTEIN APPROACHES TO 
STATISTICAL MECHANICS 

The remarkable thesis of Gibbs is that for a system in 
thermal equilibrium the observed properties of the sys- 
tem can be computed from a weighted average of the val- 
ues of the relevant observable at all possible phase points 
that lie on a constant time-slice. This approach is quite 
different from how experimental data are processed-by 
taking a time average of the 'values' of the observable at 
different times, over a very long time. The equivalence 
of this time-average to the Gibbsian ensemble average 
follows from the fascinating attribute called 'ergodicity', 
a property that is the consequence of mixing [10] . The 
ensemble average of an observable X(q,p) in equilibrium 
(indicated by the subscript 'eq' below) is defined by 

(X(q,p)) eq = Tr (p(q,p)X(q,p)) , (6) 

where 'Tr' (trace) implies an integration over the entire 
phase space in classical statistical mechanics, whereas it 
is a sum over possible eigenstates of the full 7is in Eq.(l) 
in quantum statistical mechanics. The Gibbs-Boltzmann 
weight function p(q,p) is what is called a density matrix, 
given by 

,-. exp(-/3Hs{q,p)) 

p{q,p) = y~ s ' ^> 

where f3(= (^r) -1 ) is the inverse temperature, Ub be- 
ing the Boltzmann constant. The normalization factor 
Zs , referred to as the partition function: 

Z s = Tr(exp(-pH s tf,0)), (8) 

provides the critical link between statistical mechanics 
and thermodynamics as it leads to the Helmholtz free 
energy T through the relation: 

T S = --^Z S . (9) 

From Ts all thermodynamic properties can be derived. 

It is of course outside the realm of Gibbsian statis- 
tical mechanics to address the issue of how equilibrium 
is reached. That question has to be posed in terms of 
models of noncquilibrium statistical mechanics, which are 
however not as robust and time-tested as the formula- 
tion of equilibrium statistical mechanics encapsulated by 
Eqs.(7)-(9). One model that stands out in this regard is 
based on the idea of Brownian motion [ll[. In the lat- 
ter one imagines the particle (much like the pollen par- 
ticle of Brown [12J), the Hamiltonian of which is given 
by Eq.(l), is in contact with a heat bath that drives 
stochastic (noisy) fluctuations into the system. The idea 
of Brownian motion is very physical in that if one tags the 
particle by taking camera snapshots at different times, 
its dynamics would indeed appear to be random, when 
the particle is out of equilibrium, and even when it is 
in equilibrium! The stochastic dynamics is captured 
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by the time-dependant distribution function V(q,p,t) in 
phase space that obeys the Fokker-Planck-Smoluchowski- 
Kramers equation [13| 



= {--■V, + V p .(V g V(g)+ 7 p) 



+ m~fk B TVl}V(q,p,t) , 



(10) 



where 7 is the friction constant. The quantity V plays 
the same role in non-equilibrium as p does in equilibrium. 
Thus the averaged time-evolution of the dynamical vari- 
able X(q,p) is given by 



X(t) = J dqdpX(q,p)V(q,p,t) 



(11) 



With the temperature-dependant prefactor in front of 
V 2 , it is ensured that the stationary state is indeed the 
thermal equilibrium state, described by p in Eq.(7). This 
is consistent with the fluctuation-dissipation theorem. 

Although the fluctuation-dissipation relation is a nec- 
essary condition for guaranteeing that the system tran- 
sits to the thermal equilibrium distribution, as t — > 00, 
the Brownian motion model is far from being a unique 
description for the approach to equilibrium. More signif- 
icantly, even within the Brownian motion model, there 
may be different routes to approach equilibrium. For in- 
stance, we can ask: does \m\ t ~*oo X(t) agree with (X), 
as defined in Eq.(6)? The resolution to this question 
helps our understanding of how to relate experimen- 
tally measured quantities to their theoretically calculated 
values in equilibrium, as prescribed by Eq.(6), for in- 
stance (cf., comments in the last but paragraph one in 
Sec. I). Not surprisingly then, the rich physical structure 
of the Brownian motion model has bestowed the latter 
the inspired title of the 'Einstein Approach to Statistical 
Mechanics' [Ti|. 

It is pertinent to mention here that the time-dependent 
approach, as formulated through Eq.(10), is based on 
what is called the 'Schrodinger picture'. An equivalent 
description obtains through the 'Heisenberg picture' in 
which one directly considers the dynamical equations of 
motion: 

Oq 



dt 



P 



= —muj n q 



: (qx B) 



dt -~ c 

- !P(t) + Rt) . (12) 
The set of equations (12) is called the Langevin equation 
in which the force f(t) is a stochastic noise, defined on an 
ensemble for which the distribution function is given by 
V(q,p,t). A particular realization of f(t) corresponds to 
a given trajectory, and ensemble averages are obtained by 
imposing the following constraints on the spectral prop- 
erties of f(t): 



(fit)} 







(U(t)f„(t')} = 2mjk B TS{t-t')S flv , p, v 



x,y,z . 
(13) 



III. DISSIPATIVE QUANTUM SYSTEMS 

In this section we move from the classical to the quan- 
tal domain and consider the case in which the quantum 
subsystem is put into contact with a heat bath that is also 
quantum mechanical. Before we indicate the steps nec- 
essary for Brownian motion in terms of what is referred 
to as the quantum Langevin equations Q , it is useful to 
backtrack and indicate how the classical Langevin equa- 
tions (12) themselves are derived from a system-plus-bath 
method. Here we start from a treatment of Zwanzig (l5l | 
in which the Hamiltonian in Eq. (1) is extended as 



$ , 1 2/- °jQ < 
2m j 2 J jW mjJ] 



(14) 



Upon expanding the square over the round brackets it is 
evident that the Hamiltonian contains a linear coupling 
between the coordinate q of the subsystem and the coor- 
dinate cfj of the harmonic bath with Cj being a coupling 
constant. 

From Eq. (14) it is easy to write down Hamilton's 
equations of motion, solve for the bath coordinates and 
momenta, put the solutions back in the equations of mo- 
tion for the subsystem variables and derive for the mo- 
mentum the generalized Langevin equation [l3|, [l5[ : 



mq 



-mu>Qq (<f x B) 



in I dt'q{t') 1 {t-t')+ f(t) , (15) 
'0 

where the "friction" ^(t), that appears as a memory func- 
tion, depends quadratically on Cj and the noise f(t) de- 
pends explicitly on initial coordinates and the momenta 
of the bath oscillators: 



C 2 

Rt) = E^iK-w-^?] 003 ^*) 

'—' TTljUl- 



(16) 



(17) 



Suffice it to note that Eq. (15) is exact and devoid of 
any assumption except that we have decided to integrate 
the equation of motion in the forward direction of time, 
thereby giving a sense to the 'arrow of time'. The next 
step however is a crucial one of introducing irreversibility 
by considering an initial ensemble of states, a' la Gibbs, 
in which the bath variables are drawn at random from a 
canonical distribution (Eq.(7)), yielding 



(U(t)fu(t')) = S^mkBT^t - t') 



(18) 



The final step is to go to the limit of an infinitely large 
system in order to endow the harmonic oscillator system 
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the attribute of a heat bath. Thus 



c 



, m 3 = to, Cj = —= , (19) 



where g(u>) is the 'spectral density'. Equation (16) then 
yields 

poo / \ 

7(t) = — / dw^cosfwt). (20) 

TO Jo CJ^ 

A commonly assumed form of g(co) is the one which 
yields what is called Ohmic dissipation, and is given by 



, LU < LU 



, LU > U>, 



(21) 



lu being a high-frequency cut-off. Employing Eq.(21) we 
derive Eq.(12), implying that Ohmic dissipation corre- 
sponds to constant friction 7 because the generalized fric- 
tion coefficient reduces to jS(t — t'), wherein 7 equals 



2mu 3 



If 



The discussion in the quantum case proceeds along si- 
miliar lines in which one has to however keep track of 
the fact that q and p are non-commuting operators, and 
consequently, the noise / in Eq.(17) is also a quantum 
operator Additionally, because the bath oscillators 
are to be treated quantum mechanically, the noise corre- 
lations are not 'white' , as in Eq.(13), but are character- 
ized by both a symmetric combination and a commutator 
structure, respectively given by 

({/„(*), /-(*')» = <W- r dLoM[f(LU + i0 + )]LUCOthA 
7T ,/n 2 



x cos[cj(i — t')] 

2 r°° 



m JO 
x lu sm[tu(t — t')] 



(22) 
(23) 



It is customary to rewrite Z as a functional integral [lj| : 

Z = <j)V[q 7 p, q} , pj] exp (- ^A e [q, p, q~j ,fj]j , (25) 

where h is the Planck constant and A e is the so-called 
Euclidean action, defined by 



rh0 

A e = dTC{r) , 
Jo 



(26) 



C(t) being the Lagrangian written in terms of the 'imag- 
inary time' t(= ih/3'). We illustrate in Sec. IV below the 
application of Gibbs and Einstein approaches to the cal- 
culation of the heat capacity for the charged oscillator in 
a magnetic field. 



IV. HEAT CAPACITY 

The heat capacity or the specific heat at constant vol- 
ume is the most basic thermodynamic property. It is 
defined by ^ 



(27) 



where U is the internal energy. From a statistical me- 
chanical point of view C is also related to the mean 
squared energy fluctuations given by [l8| 



C = k B (3 2 ((H 2 ) - (H) 2 ) 



(28) 



While in the Gibbs approach C can be directly computed 
from Eq.(25), employing the definition in either Eq.(27) 
or Eq.(28), the quantities (Ti 2 ) and (Tt) 2 are functions of 
the time t, in the Einstein approach. Correspondingly, C 
will also be a function of t, and the question we address 
is under what circumstances do we have the following 
equality: 



lim C(tf instcin = C Gibbs 



(29) 



At this point it is pertinent to ask: which system is (3 
(as in Eq.(22)) the inverse temperature of ? In the Ein- 
stein approach, discussed so far in this section, it is clear 
that (3 represents the harmonic oscillator bath which the 
subsytem of interest, described by Hs in Eq.(l), is ex- 
pected to eventually come to equilibrium with. However, 
because the interaction between the subsystem and the 
bath is treated exactly there is no reason for not thinking 
of the entire system, represented by the Hamiltonian Ti 
in Eq.(14), as one composite many body entity, which is 
further embedded in yet another external bath, the in- 
verse temperature of which is also given by (3\ This then 
summarizes the Gibbsian approach in which one writes 
the full partition function by replacing TLs m Eq.(8) by 
Eq.(14): 

Z = Tr(exp[-f3H]) . (24) 



A. Gibbs Approach (luq 7^ 0) 

Before we discuss the calculation of C hs for the dis- 
sipative charged oscillator in a magnetic field it is use- 
ful to indicate the steps for the simpler problem with- 
out dissipative coupling, viz; that described by Tls alone 
(Eqs.(l) and (4)) pj|. The corresponding Lagrangian for 
the two-dimensional motion in the plane normal to the 
field is given by 

11 e 
C = -m(i 2 +y 2 )--mcu 2 {x 2 +y 2 )--(xA x +yA y ) . (30) 

It is customary to work in the so-called "symmetric 
gauge" in which 

A x = ~yB, A y = ±xB. (31) 
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The Euclidean action can be written as 

A e [x,y] = y J\t[(±(t) 2 + y(r) 2 )+u; 2 (x(r) 2 + y(r) 2 ) 



- iuj c (x(T)y(T) - y(r)i(r))] , 
lo c being the "cyclotron frequency" given by 

eB 



mc 



Introducing 



x(t) = Hvj) exp(-i^r) 



(32) 



(33) 



(34) 



Turning now to the dissipative system described by the 
full many body Hamiltonian in Eq.(14) we can similiarly 
derive [l9( 

oo 4 

1 jy 



2„2 ' 
\ (43) 

where 7(1^) is the frequency (ie., Vj) — dependent friction 
coefficient. The Ohmic dissipation model, discussed ear- 
lier in Eq.(21) that yields constant friction, is not suitable 
for calculating Z as it leads to a singularity. In order to 
regularize the latter it is convenient to introduce a 'Drude 
cut-off' by writing the spectral density as (cf., eq.(21)) 



where Vj's are the so called Matsubara frequencies, de- 
fined by 



2^j 

h(3 



j = 0,±l,±2, 



we find 

A e [z+,Z- 

where 



+ (v 2 +loI - iuj c i/j)z^(Pj)z^(Vj)] , (36) 



z±{vj) = ± iy(vj)) 



(37) 



As shown in Ref. [19] the partition function Zg in equa- 
tion (8) can be written as (cf., also Eq.(25)) 



.7=1 



(38) 



where, 



Z+ = 



1 r 00 

\ dz+(Q)exp 
7rnrp/m J -«> 



00 />oo />oo 



dR,ez+dhnz + 



00 j — 00 



exp [— m(3 (i/j + 



and 



Uln — UjJ c V 



) (Rez 2 + + lmz 2 + )] , 
(39) 

(40) 



Carrying out the Gaussian integrals we find 
_ 1 



Hence, 



/ 1 



f3faj a (v 2 + uj' 2 - iu c Vj) 



n 



j=l (v* + wg) + u*u. 



2, ,2 

3 



(41) 



(42) 



2m7 



7TC 1 



(35) 

Correspondingly (cf., Eq.(21)), 



(44) 



S • (45) 



All our results in the sequel are restricted to Ohmic- 
Drude spectral density (Eqs.(21) and (44)), though it is 
known that other forms of frequency-dependence of the 
spectral density yield diverse forms of power-law depen- 
dence of the specific heat at low-temperatures [20j | . 

Inserting this form of the friction coefficient in Eq.(43) 
the internal energy U can be calculated as 



p v v 



V V V V 



(46) 



where ip(z) is the digamma function and the arguments 
are: 



Ai + A2 + A 3 = ojd + i<jJ c , 

A1A2A3 = WqW£) 



AiA 2 + A 2 A 3 + A3A1 = uj 2 + ryuj D + ilo c uj d 



(47) 



The corresponding primed A's are obtained from the com- 
plex conjugate of Eq.(47). Finally, it is easy to derive for 
the heat capacity the expression (cf., Eq.(27))[l!| 



Ogft = -2** + **£{(£)V(£) 



v ] v 



(48) 



We are now ready to discuss the low and high- 
temperature limits of the heat capacity. 
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(a) Low- T limit 



where 



i^fGibbs 



2tt 7 k%T 



8tt j 



fr 4 

" B 



15 (S-wo) 3 



af T 3 + C(T b 



3(w*+wg) 



U!q L0 B LUq L0 B ' 



(49) 



the partition function for a charge in a magnetic field 
(without the oscillator potential). In this case only two 
roots Ai and A2 (cf., Eqs.(46)) matter [l9| and we find 



Z{ujo = 0) 



Nmfi 2 2 

(7 +^c>- 



8tt 3 

The heat capacity becomes 



y~*Gibbs 



k B +k B y^A 

k=l 



(7(^)) 



v v 



(52) 



Curiously, to leading order, the presence of the mag- 
netic field through the cyclotron frequency disappears 
from C^^o) ' ^ ne expression of which matches with that 
of a two-dimensional quantum oscillator (Einstein oscil- 
lator). The result in Eq.(49) has been much in discussion 
in recent times, in the context of the third law of thermo- 
dynamics as it provides a satisfactory power-law behavior 
in temperature [4|. 

(b) Hight -T limit 

At high temperatures (huj Cl Huiq, try, ftw B << k B T) 
our quantum system is expected to be described by clas- 
sical statistical mechanics. We find 



-.Gibbs 



J^2 



(50) 



where 



12k, 



-{cu 2 + 2Lu 2 + 2~fu D ) 



In the limit of infinite temperature, therefore, we re- 
cover the expected 'equipartition' result: 



(51) 



where the factor of 2 comes from two dimensions, each 
of which contributes k B to the specific heat, \k B arising 
from the kinetic energy while the other half from the 
potential energy. 



B. Gibbs Approach (ujc, = 0) 

While studying dissipative Landau diamagnetism we 
have learnt that taking loq = at the outset yields puz- 
zlingly different result from keeping loq fixed, evaluat- 
ing the partition function, calculating its derivatives and 
then setting loq — @ . It is already evident from the low- 
temperature specific heat (Eq.(49)) that it is not mean- 
ingful to take the limit of ujq = without 'fixing' the 
coupling with the heat bath characterized by the friction 
coefficient 7 ! It is therefore of interest to take a relook 
at the heat capacity calculation by investigating afresh 



V j V 



2k, 



\ V ) 



(53) 



We now discuss the low and high temperature limits of 
Eq.(53). 

(a) Low- T limit 

Using asymptotic expansions as before, we find 



27r 7 (l-^M 2 



y^fGibbs 



3 h^+ul) 

4 



kiT 



(af - af)T 3 + 0(T 5 ) . (54) 



where 



8tt 3 



i.i 



( 7 3 - 3 7 ^ c 2 ) _ 37 



15 ft 3 V(7 2 + ^ c 2 ) 3 I V(7 2 + ^ c 2 ) 3 
(w 3 - 3uj c j 2 ) f u c 3 



, (^) 3 + 3 7 ^- 

V(7 2 +^ c 2 ) 3 V W ^1 



a _ 8^ 3 k% 



4 15 {hw D f 

While Eq.(54) is in conformity with the third law of 
thermodynamics with identical linear temperature de- 
pendence as in the case of loq 7^ 0, but, is free from the 
singularity issue in Eq.(49) (for u>o = 0). It leads, in the 
limit of u>d = 00 (infinite Drudc cut-off) to the result: 



7 



^tGibbs _ ^Tr 2 

Further, for very strong magnetic fields (7 << cj c ), 

^Gibbs _ 2?r 7 k 2 B T 
°^=°) ~ 3 h ' 



(55) 



(56) 



a harmonic oscillator like result with the cyclotron fre- 
quency uj c replacing luq. On the other hand, for weak 
magnetic fields (7 >> uj c ), 



y-<Gibbs _ 2-7T k B T 1 
°^=°) "3 h 7 ' 



(57) 
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the free particle result in which the friction coefficient 7 
appears in the denominator, in agreement to the corre- 
sponding result given in 21], after a proper counting of 
the degree of freedom. 



(b) Hight -T limit 
We find 



l2k B T 2 



{lo 2 c + 2^0) 



(58) 



Again, equipartition theorem for a free particle (in 2 di- 
mensions) prevails at T = 00. 

Thus the classical limit of the Landau problem, as far 
as the heat capacity is concerned, is that of free particle 
whereas an additional (parabolic) constraining potential 
yields harmonic oscillator behavior. 



C. Einstein approach (ujo 7^ 0) 

We will now focus on the Einstein approach based on 
the Langevin equation (15)which can be recast into the 
following convenient form [9(: 



2„ F (t) 



z + / dt'j(t - t')z(t') + Jiz = 
Jo m 



(59) 



where 

z = x + iy, F = f x + if y , and 7(f) = 7(t) + ilu c . (60) 

In order to find the time-dependent specific heat we need 
the internal energy which is the statistical average of the 
Hamiltonian given by 



Ti. = —razz* Huj c H — muxkzz^ . 



(61) 



We therefore need the equal-time correlation functions: 

G(i) = (*(*)**(*)) , (62a) 
C 2 (i) - (z(t)zHt)) • (62b) 

The correlation functions in Eq.(62) can be found from 
the analytic continuation to t' = t of the unequal time 
correlation functions, eg., 



d(t,t') = (z(t)zHt')) , 



(63) 



where z(t) can be further expressed in terms of the re- 
sponse function x(i) as 



z(t) 



I dr X (t-r) 
Jo 



F(t) 



(64) 



The former is the inverse Fourier transform of x( w ) that 
can be easily written from Eq.(59) as 



1 



2-7T (— LU 2 — ILOJ + LOq) ' 



(65) 



with 



7(w) = iuj c + 7(w) = ilu c + 7 — — . (66) 



From Eq.(63), 



Ci(M') = / dr [ dr'x(t-T) X *(t' -r') 
Jo Jo 

„ (F(t)FHt')) 



where 



m 



+ OC 



(F(r)Ft(r')) = / ^/(^e^— ') 



(67) 



(68) 



with 



m = - 7 T?~ 2 ^ w^) - 1] ■ (69) 



tlLO 



7T (UJ 2 D + UJ 2 ) 



"2kT' 



Our strategy is to first calculate £i(t, i') and C,2(t,t') 
(for details, see the Appendix A), then set t = t' and fi- 
nally, in order to extract the thermal equilibrium internal 
energy E, take the limit t = 00. We find 

E = (H) = - \tiw c + \mu>l lim Ci(*) + \ m Hm Ca(*) 

Z Z t — >oo Z t — >oo 

h ^ ( A 
= 2k B T + — U(l + -f )[2wfo + Pj ] 

3=1 



+ v(i + ^)[2w 2 <z;+ P ;" 

where 



Pj 



Xj[jojp - iuj c (Xj - uj d )] 



(70) 

(71a) 
(71b) 



In the denominators of Eqs.(71), the notation FT., implies 
that the j = j 1 terms are excluded from the product. The 
quantities q'j and p'j are obtained by priming the A's, the 
latter having been already defined in Eq.(47). 
Finally, the equilibrium specific heat is given by 



Einstein 



dE 
df 



Ztt * — ' I V V 

3=1 



V V 



(72) 



where tp'(z) are the trigamma functions fl9j ]. 

We may now discuss the low and the high temperature 
limits of Eq.(72). 
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(a) Low-T limit 

Employing the asymptotic expansion of the digamma 
function: 



^ (Z) ~ z + 2z 2 + 6z 3 30z 5 



(73) 



we find 



C 



Einstein 

LU^O 



2tt 7 k 2 B T E 



3 ujI h 



+ af T 3 + (9(T 5 ) . (74) 



where 



a, = 



8^ 3 7 4 f3K 2 +^ 2 ) 
15 w (M)) 3 



27 W 



w 



Wfl 7^0 ^0 



As required by the third law of thermodynamics the spe- 
cific heat does vanish as a power law as T — > 0, exactly 
in the same manner as in the corresponding Gibbs ex- 
pression (cf., Eq.(49)), but interestingly the coefficient of 
the next higher order term (oc T 3 ) differs from the Gibbs 
result. 



As discussed in Ref.[19], one of the three roots, viz. Ai 
vanishes for lo = 0. Consequently (see Appendix B, for 
details), 



E(uj = 0) 



k B T 

h_ 

2^ 



te^(l + ^)+^(l + ^) 



+ p' 2 ^(l + -)+P3^ + -)} ■ (79) 



As before, the derivative of E with respect to tempera- 
ture yields an expression for the specific heat in terms of 
the digamma functions, which can be further analyzed in 
the low- and high-temperature limits. 



(a) Low-T limit 

Again, using the asymptotic expansion of the digamma 
function (cf., Eq.(73)), we find 

c Einstein = I fc| T - af T 3 + 0(T 5 ) . (80) 



where 



3 H^ + ujI 



(b)High-T limit 



At high temperatures, 



where 



CEinstein 



a 2 



(75) 



8tt 3 



a q = 



Kg 



+ 10( 



(7 3 - 3 7 ^ 2 ) 

15 ftV(7 2 + ^ c 2 ) 3 1 v% 2 + 

Ue ,2 7(7 2 + W 2 ) 



12fc s 



(w 2 + 2^0 +7w D ) 



At infinite temperatures the classical equipartion result 
is restored. But again, in the next higher order term (in 
-ly), the Einstein result differs from the Gibbs result by 
a cut-off-dependent term: 



r Einstein 



C Gibbs 



l2k B T 2 ' 



(76) 



While the expression in Eq.(80) is in conformity with the 
third law of thermodynamics, as expected, it differs from 
the corresponding Gibbsian result of Eq.(54) in terms of 
different dependencies on the Drude cut-off wd! Apart 
from this issue the strong and weak magnetic field cases 
follow the behavior discussed earlier, below Eq.(54). 

(b)High-T limit 



D. Einstein approach (u>o = 0) 

We now return to discuss the Einstein result for the 
specific heat due to the presence of the magnetic field 
alone, ie., in the absence of the parabolic well. The rele- 
vant Hamiltonian is 



and hence 



E 



-—hw c 
2 



-m lim [C 2 (*)U=o 

Z t^oo 



(77) 



(78) 



CEinstein ;„ 
uj =o — Kb - 



n 2 



12k B T 2 



(uj 2 c + ju D ) . (81) 



Finally, in the high-temperature limit, equipartion result 
obtains, but once again, there is a correction term over 
and above the Gibbs result that is cut-off dependent, as 
we found earlier in the wq 7^ case in Eq.(76): 



j-2 

s-tEinstein /-iGibbs , 7^-P 



(82) 



where C2„=q is given by the high-T expression in Eq.(58). 
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w c # 0, wo / 


w c / 0, u; = 


U c = 0, WO = 




Low 

Temperature 


High 

Temperature 


Low 

Temperature 


High 

Temperature 


Low 

Temperature 


High 

Temperature 


Gibbs 
Approach 


2vr 7 k ' 2 B T 
3 h 

-af T 3 + 0(T 5 ) 


2fc s -^ 


(1-—) o 
1 U / » ^ 1 




2tt fc B T q _ 7 \ 
3 ft 7 ujj ' 

-(a|fU=o -af)T 3 
+0(T 5 ) 


/ . . ^ V^U — u j — ■ / 


3 ft 7 2+^ KbJ 

-( a3 G -a 4 G )T 3 
+0(T 5 ) 




Kb t 2 


Einstein 
Approach 


2vr 7 fc|T 
3 ft 

-af T 3 + 0(T 5 ) 




2ir 7 1 j„2 




+0(T 5 ) 


y. £ af |(w c =0, w =0) 


3 B^KflJ 

-af T 3 + 0(T 5 ) 







TABLE I: Comparison of Specific Heat in the Gibbs Approach and the Einstein Approach in different limits. 





Specific Heat 


Magnetization 




Low Temperature 


High Temperature 


(7-0) 


LOa — > 0, t — > 00 


^i^klT-aiT" + 0(T 5 ) 




-Scoth(^) 


t — ► oo,a;o — > 


Singularity 


2fc s a ^r° 


^[^Z-coth(^)] 



TABLE II: Specific Heat and Magnetization in the limit of vanishing confinement frequency in two sequences. 



V. SUMMARY 

Summarising, we study the various limiting behavior of 
the specific heat of a dissipative charged harmonic oscil- 
lator in a uniform magnetic field, obtained from the par- 
tition function approach (Gibbs' method) and from the 
steady state of corresponding quantum Langevin equa- 
tion (Einstein's approach). The specific heat obtained 
from both these methods shows linear T dependence at 
low temperatures, which is in agreement with the third 
law of thermodynamics. At high temperatures the spe- 
cific heat approaches a constant value depending on the 
number of degrees of freedom of the system. Although, 
both the Gibbs and Einstein approaches are in confor- 
mity with the third law of thermodynamics and the 
equipartiton theorem, at low and high temperatures re- 

I 



In Table I, we summarise our results for the Specific Heat 
in different limits for both the Gibbs and Einstein ap- 
proaches. In the limit of lojj — > 00, both the Gibbs and 
Einstein approaches give the same thermodynamic re- 
sults. However, for a finite cutoff frequency ujd, the re- 
sults differ in next to the leading order at both high and 
low temperatures. The results summarized in Table I 
lead to the following conclusions : 

I. At low temperatures the specific heat is linear in 
temperature and hence the dissipative environment 



spectively, they differ from each other in detail, beyond 
the leading order. In the limit of vanishing confinement 
frequency (uio — > 0), the specific heat of the oscillator 
becomes singular at low-temperatures and manifests ex- 
tra degrees of freedom counting at high temperatures. 
The specific heat of the free particle cannot be obtained 
from the equilibrium value (t — > 00) of the specific heat 
of the oscillator just by taking the ojo — > limit. It is 
evident that the order in which one takes the t = 00 and 
ujo = limits yield qualitatively different answers for the 
specific heat. While in the Einstein approach, the free 
particle-like specific heat emerges by taking the u>o = 
limit first before considering the t = 00 limit, the Gibbs 
approach is plagued by a singularity issue, for uio = 0, in 
the low-temperature limit (cf., Eq.(49)). 



I 

restores the third law of thermodynamics. 

2. In the presence of the oscillator potential, the low 
temperature behavior of the specific heat goes as 

and is therefore singular in the limit of u>q — > 
0. Thus the results of the unconfined particle can- 
not be recovered in this limit. 

3. The high temperature specific heat approaches a 
constant value independent of the confinement po- 
tential and depends only on the number of degrees 



10 



of freedom in agreement with the equipartition law. 
Again, the results of the unconfined system cannot 
be recovered in the limit of vanishing confinement 
frequency No- 



where XjS and X'jS are given by the Vieta equations 
(Eq.(47)). We can write Eq.(67) as 

i r + °° 

Ci(t,0 = (*')) = / dam 



While the issue of recovering the results of the unconfined 
particle, starting from the confined system and taking the 
limit of vanishing confinement frequency loq cannot be 
resolved at the equilibrium level, the Einstein approach 
has the intrinsic advantage of obtaining the results in the 
process of equilibration. The equilibrium results can be 
arrived at by taking the limit of t — > oo. Hence, one 
could in principle ask the question, what would happen 
if the confinement frequency uio is taken to zero, before 
the limit t — > oo is taken. A similiar result was ob- 
tained for the case of a particle in a harmonic oscillator 
potentially. The results for the two different sequences 
of taking the limits is summarised in Table [ill It is clear 
from the table that, if the limit of uio — ► is taken be- 
fore the limit of t — > oo, one can actually recover the 
results of the unconfined system for the specific heat and 
magnetization. It is curious to note that the result for 
magnetization obtained from this sequence of taking the 
limits is inconsistent with the Landau results, whereas 
when the limits are taken in the other way round, the 
Landau result is recovered. This is, however, due to the 
fact that the Landau result for magnetization can only 
be recovered in the presence of a confinement potential. 
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• / dux{u)- — " 



hoo f„iuit' „iu>'t'\ 

x / dwYV)- 



^ (e m -e iuit ) 
-i(u>' — ui) 



The two integrals, defined by 



h = 



^g — iCut ^—iu)t\ 

i(uj — Cj) 



(A3) 



(A4) 



duW V ■/ i ~\ > (A5) 

can be expressed as 

j = 2tt (A! - ui D ){X 2 - A 3 )(e~ iat - e-^') 
1 iA^ (w + iAi) 

(X 2 -LO D )(X 3 -X 1 )(e-^ t -e-^ t ) 



(u> + iX 2 ) 
(X 3 -m D )(X 1 ~X 2 )(e- lM ~e-^ t ) 
(u> + iX 3 ) 



} , (A6) 



= 2n (X[-Lo D )(X! 2 -Xi a )(e^' -e-^') 
2 iA ,{ (w-tAi) 



(X' 2 -m D )(X' 3 -X' 1 )(e^ t -e~^) 

{Co - iX' 2 ) 
(A^-^XA'i-A^Xe^'-e-^') 

(Co - iX' 3 ) 



} • (A7) 



where 
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useful discussions. SD is grateful to the J. C. Bose Fel- 
lowship of the Department of Science and Technology for 
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A = (Ai — A 2 )(Ar - A 3 )(A 2 - A 3 ) , (A8) 
A' = (A i -A 2 )(A' 1 -A 3 )(A 2 -A 3 ) . (A9) 



Eq.(A3) then yields 



APPENDIX A: EINSTEIN APPROACH (w + 0) 



With the help of the Drude cut-off frequency we can 
write x( w ) as 



X(w) = 



(ojd — iio) 



[iLU 3 — L0 2 {iOD + iiO c ) — ILU^LUD + iui c UlD + Wq) + IOqLJd] 

(Al) 

Alternatively, 



xM 



{ui + iui D ) 



(u + iX±)(uj + iX 2 )(u> + iX 3 ) 



(A2) 



1 



4ir 2 m 2 
1 

47r 2 m 2 



dCof(Co)hI 2 



— oo 

+oo 



dCo — - — 



7r (uj 2 d + Cj 2 ) ' 



x coth (^) /l/2 



1 



47r 2 m 2 



dCo — ; J c . n hOIi I 2 



7r (ujjj + Cj 2 ) 



(A10) 



The second integral vanishes for symmetry reasons, so 
that only the integral containing cotangent hyperbolic 
contributes. In order to find out the equal time correla- 
tion function we set t = t' . In that case the coeffi- 
cients of e~ ZUJt and e lut matter, because in the product, 
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these are the only time independent parts. Now substi- 
tuting Eqs.(A6) and (A7) in (A10), we can easily sepa- 
rate the mean squared average into two parts, one that 
is completely time independent and the other which is 
an exponentially decaying (time dependent) one. In the 
limit of t — > oo, the time dependent parts vanish and we 
are left with the equilibrium value. Finally, 



&(*) = 



47T 3 TO 



/+oo 
doo 
-oo 



x coth(^)^ 



where 



We can write 



2tt 



(lod — iu))e 



-iut 



i (Co + + i\2)(u> + 1X3) 

2tt (lo d + iCo)e lM 

i (Co — iA^)(w — i\' 2 ){ti — iX' 3 ) 



(All) 



(A12) 
(A13) 



We will calculate &(t, t'), which is defined as 



fc r+00 

Ca(*, - (z(t)z\t)) = / dwwY'cothf 

run J_ oc 

J —OO 



run 

2kT Hlo, 



fko 
2kT> 

(A17) 
A' 



m 



7 = 1 



- ^{E^E — x\ 

*• 7=1 n=l 71 + 17 7=1 n=l n + -f- J 



3 oo 



+ 



(A18) 



where qj and the ^ are defined in Eq.(71a), and pj and 
p'j are given by (71b). We now use a transformation Pj = 

Pj + i^, such that Y^j=i Pj = 0' since X^j=i Pj ~ ~* w c- 
Therefore 



Ci(t) = <*(*)**(*)> =Qi-Qa 



where 



Qi 



Q2 = 



2nm 



/+00 
-00 



dojcothl— — ) 
y 2kT J 



(lo d - iCo) 



(Co + iAi)(w + iA 2 )(u) + iA 3 ) 



ft 



27rm 



/+00 
-00 



, , hoo . 

aucoth — — 
v 2fcT y 



(cj — iAi)(w — i\' 2 )(Co ~ ^3) 



(A14) 



(A15) 



Assuming that the time is long enough compared to the 
relaxation time, we can ignore the integrals containing 
Ai, A2, A3. After simplifications 



d(t) = (z(t)zHt)) 



2kT 



mcon 



+ 



^j^{ qji) (i + h ) + q ^( l + %\ 

mn ^— f I zv J z/ I 



(A16) 



where ^(1 + Zj) is a digamma function, ^ = 27T ^ T , and 
the and the <^ are defined in Eq.(71a). We can ob- 
serve from Eq.(A16), (since (r 2 ) = (z(t)z*(t))) that the 
equipartition theorem is satisfied for this two-dimensional 
problem. 



3 \ J 3 - 1 

n H — - « J — 3 - 



3 



\ .L0 C \ ^ 1 



3 00 1 

j=l n=l 72 + ;y 



(A19) 



In a similiar fashion wc can use a transformation P- 

p'j-i^, in such a way that X^=i = since Z)j=i 
2w c , hence 



3 00 



3 00 



e^e-^ = e^e-^ 

j=1 «=Hl+f i=l n=l n + 
3 00 1 

+ E*yE^ 

j=i „=m + 

= + 

3 00 1 

+ E*tE^t-(A2o) 



n=i n + -f 



Substituting Eqs.(A19) and (A20) in Eq.(A18) and using 
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three important properties of the digamma functions [2c 



find 



ip(x) - ip(y) = 



- y) 

xy 



oo 1 

£[— 



^— ' n + v n + x 

n—l ,J 



+ z) = lj)(z) + - , 



N 



a,- — - — 

f-f ^ 71 + 2j 

j— 1 n—l J 



N N 

(A21) 



we obtain 



C2(*) 



2fcT 



r/v 



(i(t)it(t)) 

^E{^(i + ^) + ^(i + ^)j 

i=i k J 
run ' ^ 7 z/ 

5=1 k J 



ftw c 
m 



(A22) 



From Eq.(A22), we can calculate the mean squared av- 
erage of the kinematic momentum of the particle in a 
magnetic field, given by 



Einste 



2k f 



.7 = 1 k 



+ g ;^v/(i + ^) 



.7=1 k 



(A25) 

where are the trigamma functions and kg is the 

Boltzmann constant. Finally employing the recurrence 
formula for trigamma functions leads to 



4>'{l + z) = ip'(z) 



1 



and also 



3=1 
3 



l 



3=1 A ' A ' 

from which we obtain Eq.(72). 



(A26) 



((P--i*) 2 ) = m 2 (z(t)zUt)) -mhu. 
c 

= 2mkT 

,2 3 



mhuj7{ \ -\ ( , . A,- , , . A 

+ — - J2 1^ 1 + — ) + + - 

3=1 k 

771 ft f , , A.,, , , . A' 



(A23) 



In the limit of a vanishing magnetic field, the two aver- 
age values which we calculate are similiar to the result 
obtained for a damped harmonic oscillator, as given by 
Weiss Q, of course with a different degree of freedom. 

The internal energy can be obtained as 



E(uj ) = (H) = ^m(iit) - l -hw c + ^mLul(zz^) . (A24) 



Taking the derivative with respect to temperature, we 



APPENDIX B: EINSTEIN APPROACH (w = 0) 



Here we provide details of the calculations for the case 
of luq = 0. Here, one of the three roots, viz., Ai vanishes 
and we are left with just two roots. From the Vieta equa- 
tions given in Eq.(47), we can write the new equations 
for this particular case as A2 + A3 = ujd + i^c , A2A3 = 
w/j(7 + iw c ). In the limit of vanishing harmonic oscillator 
frequency, the energy is obtained as Eq.(78) 



E = -\hw c + ~m lim [C 2 (i)k=o • (Bl) 



We can write £ 2 (i) 



zz>) as 



hm C2W wo=o = H 

t^oo m m 



7717T y V V 

+ + -(B2) 
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where 



[yu> D - iLu c (X 2 -u) D )] 
P2 = — and 



(A 2 - A 3 ) 
Yyujp - iuj c (X 3 - m D )] 
(A 2 - A 3 ) 



(B3) 



The primed roots are calculated from complex conju- 
gates. Hence, the internal energy is 

E(lu q = 0) = k B T 



+ ^L(l|^)l|.3# + ^) 

2tt { v v 

+ P 2 V>(l + ^)+P 3 V>(l + ^)} • (B4) 



Correspondingly, the specific heat becomes 



z-iEiristein 



= -k B 

2n I v v v v 



This form of the specific heat has been used in the text 
as the basis of our discussions of the low and high tem- 
perature limits, via Eqs.(80) and (81). 
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